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We study a granular gas heated by a stochastic thermostat in the dilute limit. Starting from 
the kinetic equations governing the evolution of the correlation functions, a Boltzmann-Langevin 
equation is constructed. The spectrum of the corresponding linearized Boltzmann-Fokker-Planck 
operator is analyzed, and the equation for the fluctuating transverse velocity is derived in the 
hydrodynamic limit. The noise term (Langevin force) is thus known microscopically and contains 
two terms: one coming from the thermostat and the other from the fluctuating pressure tensor. At 
variance with the free cooling situation, the noise is found to be white and its amplitude is evaluated. 

PACS numbers: 



I. INTRODUCTION 



Typically, a granular system is defined as an ensemble of macroscopic particles which collide inelastically, i.e. 
part of the kinetic energy of the grains is dissipated in a collision. This simple ingredient gives rise to a very rich 
phenomenology which is of interest not only from practical or industrial perspective, but also because of the resulting 
new theoretical challenges [J 0, H, Q . One of the most widely employed idealized model for granular fluids is a system 
of smooth hard spheres (or disks in two dimensions) whose collisions are characterized by a constant coefficient of 
normal restitution 0, [||. For this model, and considering that the particles move freely between collisions, kinetic 
equations have been derived: starting from the dynamics of the particles, it is possible to derive the corresponding 
Liouville equation, and the Boltzmann equation results in the low density limit [7|, |8[ . This kinetic equation has been 
extensively used to address many fundamental questions such as the derivation of the hydrodynamic equations, with 
explicit expressions for the transport coefficients, which have been derived by the Chapman-Enskog method 0, IToj 
and also via the linearized Boltzmann equation Due to the inelasticity of the collisions, the total energy of an 
isolated granular system decays monotonically in time. In the fast-flow regime, it has been shown numerically that, 
for a wide class of initial conditions, the system reaches the so-called Homogeneous Cooling State (HCS), in which all 
the time dependence of the one-particle distribution function goes through the granular temperature, which is defined 
as the second velocity moment of the distribution [l2|, [H[ ■ This state has been extensively studied in the literature 
and very recently the fluctuations of the transverse velocity have been analyzed [3, [HI. It has been found that 
the transverse velocity fulfills a Langevin equation but, in contrast to the elastic case, the noise is not white and the 
second moment of the fluctuations is not only controlled by the viscosity but also depends on a new coefficient. Similar 
results are found for the other hydrodynamic equations [1 61 ] - The study of fluctuations in the HCS is important for 
the development of a general theory of fluctuations in granular systems because it defines the reference state from 
which macroscopic hydrodynamic equations can be derived @. In this sense, the HCS plays, for inelastic gases, a role 
similar to the equilibrium state for molecular gases. 

On the other hand, there are situations in which the grains cannot be considered to move freely between collisions. 
If, for example, the grains are immersed in a medium which acts as a thermostat, the system may reach a stationary 
state in which the energy injected by the thermostat is compensated by the energy dissipated in collisions. Note 
that, if the grains are Brownian particles, the interstitial medium injects energy into the granular system, but also 
acts as an energy sink due to frictional forces. One of the simplest mechanism that can be considered to thermalize 
the system is a white noise force acting on each grain, which results in the so-called stochastic thermostat [13, HH, 
El lEjp, [H, HI HI HE HI HI HI. One important point is that the distribution function differs from that of the 
HCS [131, an d that it is this distribution which plays the role of the "reference state". The non-equilibrium steady 
state that the system reaches in the long time limit exhibits long-range correlations which are in agreement with 
the predictions of fluctuating hydrodynamics [20| . The latter description was introduced phenomenologically, and is 
expected to be valid in the vicinity of the elastic limit only. The objective of this work is to derive these equations 
from a more fundamental point of view and without the restriction of small inelasticity. More precisely, we adapt the 
formalism worked out in 15| for the free cooling, to the present driven case. Starting from a Boltzmann-Langevin 
description, we derive a fluctuating equation for the transverse velocity identifying the noise of this equation. Under 
certain hypothesis to be clarified in the text, we obtain that the correlation function of the noise is well approximated 
by the one introduced in Ref. [20| , where the internal noise contribution (excluding the "external" noise term directly 
stemming from the thermostat) fulfilled a fluctuation-dissipation relation as for conservative fluids [29l ]. 
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The remainder of the paper is organized as follows. In Sec. [II] previous results for a system heated by a stochastic 
thermostat are presented, such as the equations for one-particle distribution function and the two-particle correlation 
function. In Sec. IIIII the Boltzmann-Langevin equation for this system is derived and the properties of the noise are 
inferred. The particular case of the transverse velocity field is analyzed in Sec. IIVI and finally, the conclusions are 
presented in Sec. [V] 



II. STOCHASTIC THERMOSTAT: PRELIMINARY RESULTS 



The system considered is a dilute gas of N smooth inelastic hard particles of mass m and diameter a. The position 
and velocity of the ith particle at time t will be denoted by Ri(i) and Vj(t), respectively. The effect of a collision 
between two particles i and j is to instantaneously modify their velocities according to the collision rule 

V'j = v i + ±±^(*.v y )*, (1) 

where V,j = Vj — ~Vj is the relative velocity, <r is the unit vector pointing from the center of particle j to the center 
of particle i at contact, and a is the coefficient of normal restitution. It is defined in the interval < a < 1 and it 
will be considered here as a constant, independent of the relative velocity. Between collisions, the system is heated 
uniformly by adding a random velocity to the velocity of each particle at equal times. The driving is implemented in 
such a way that the time between random kicks is small compared to the mean free time. Then, between collisions, 
the velocities of the particles undergo a large number of kicks due to the thermostat. In addition, we will assume that 
the "jump moments" of the velocities of the particles verify 

Bij ^ = Hm (A ^ A A ^ ) "°- = e^5 01 + f - 1)6^ (2) 
i,j = l,...,N and /?, 7 = 1, . . . , d 

where we have introduced AV^ i( g = Vi^it + At) — Vi^it), V^p(t) being the /? component of the velocity of particle i 
at time t. We have also introduced the strength of the noise, £q, and {. . . ) n oise, which denotes average over different 
realizations of the noise. The non-diagonal terms (corresponding to i ^= j and /3 = 7) are necessary in order to 
conserve the total momentum [321 . 



A. Kinetic Equations 



Given a trajectory of the system, one-point and two-point microscopic densities in phase space at time t arc defined 

by 



N 



F l (x 1 ,t) = ^2S[x l -X i (t)], (3) 

i=l 

and 

N N 

F 2 ( Xl ,x 2 ,t) =J2J2 6 ^ Xl ~ MtMxi - Xjit)}, (4) 

i=l j=£i 

respectively. Here Xi(t) = {Ri(t), Vj(t)}, while the Xi = {r^, v^} are field variable referring to the one-particle phase 
space. The average of F\(x\,t) and Fz (xi , x% , t) over the initial probability distribution of the system p(T,0), where 
r = {X\(t), . . . , Xpjit)}, are the usual one-particle and two-particle distribution functions, 

/i(as 1 ,t) = (-Fi(a:i,t)>, Mx!,x 2t t) = (F 2 ( Xl ,x 2 ,t)), (5) 

with the notation 

(G)= [ dTG(T)p(T,0). (6) 
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In the dilute limit, assuming molecular chaos, i.e. that no correlations exist between colliding particles, and that 
the sizes of the jumps due to the thermostat are small compared to the scale on which the distribution varies, the 
equation for the one-particle distribution function, fx (x% ,t), is the Boltzmann-Fokker-Planck equation [13] 



Cl C\ (-1 £)2 

-/!(x!,t) + Vl • g^flM = J[h\h] + f Q^fl(xi,t), (7) 



where 

= y^2^(ri 2 )To(v 1 ,v 2 )/ 1 ( a ; 1 ,i)/ 1 (x- 2 ,i) (8) 

and 

T (v u v 2 ) = a*- 1 J d&0(& -v 12 )(& - l], (9) 

is the binary collision operator. The operator &T 1 changes the velocities to its right into the pre-collisional velocities 

v* = vi - -^(ct- Vi 2 )ct, (10) 

„ 1 + a . A . A 

V 2 = v 2+-^ (<T-Vi2)<T. (11) 

As can be seen, the last term in does not appear in the free cooling case and depends on the strength of the 
heating, £ . 

Let us introduce the two-particle correlation function through the usual cluster expansion 

f 2 (xi,x 2l t) = f 1 (x 1 ,t)f 1 (x 2 ,t) + g 2 (x!,x 2 ,t). (12) 
Neglecting three-body correlations, the equation for the correlation function g 2 (xi, x 2 , t) was derived in [32| and reads 



d d d 

777 + Vi • h V 2 • ^~ 

at or i or 2 



g 2 (xi,x 2 ,t) = 5(ri 2 )f (-vi,v 2 )fi(xi,t)fi(x 2 ,t) 

({x 2 ,t)]g 2 (x 1 ,x 2 ,t 
fi(xi,t)f 1 (x 2 ,t), 



+ [K(x u t) + K{x 2 ,t)]g 2 {x 1 ,x 2 ,t) 



N dv\ dv 2 ' 



(13) 



where 



K( Xi ,t) = J dx 3 5(r i3 )fo(v i ,v 3 )(l+V i3 )f 1 (x 3 ,t) + ^-^, (14) 
with V a b an operator that interchanges the label a and b in the quantities to its right. 

B. The stationary state 

It has been shown numerically that, for a wide class of initial conditions, the system reaches a homogeneous 
stationary state [20| in which the energy input from the thermostat is compensated by the energy lost in collisions. 
In this case the Boltzmann-Fokker-Planck equation reads 

f £sfH(y l ) + J\f II \f a ] = 0. (15) 

For the sake of simplicity, let us introduce the following dimensionless distribution 

Mvi) = ^Xh(c), (16) 
v d H 
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with Xh(c) an isotropic function of the modulus of c, vh = (^ L ) 1 ^ 2 5 c = T~ anc ^ Th the granular temperature 
defined as 



d f 1 

-n H T H = / dv-mv 2 f H (v). 

By introducing equation (fll)|) in equation ()15[) . we obtain a closed equation for \h 



I 2 9 

dc 2 T(ci,c 2 )xii-(ci)xij(c2) + y^2^( c i) = °' 



where 



r(ci,ca) = / do-0(ff •cia)(o--cia)[o- 2 6l 1 - 1], 



(17) 



(18) 



(19) 



is the dimensionless binary collision operator and f 2 = %^ is the dimensionless strength of the noise with i = 

(njji7 d_1 ) _1 proportional to the mean free path. In the case of the correlation function, it is convenient to introduce 
its dimensionless counterpart, gn , as 



92,h(xi,x 2 ) 



^jffif(ll2,Ci,C 2 ), 



Pv 



(20) 



u 



where we have introduced the dimensionless length scale I = r/£ and I12 = li 1 2 . The reduced distribution fulfills 



d 

A( Cl ) + A(c 2 ) - c 12 ■ — 
CI12 



<?ij(ll2,ci,c 2 ) = S(1 12 )T(c 1i c 2 )xh{ci)xh(c2) 
~ 2 n H i d d d 



N dci dc 2 



Xh{ci)xh{c 2 ), 



where A(cj) is the linearized Boltzmann-Fokker-Planck operator 



£2 Q2 



A(ci)h(ci) = / dc 3 T(c ? ,c 3 )(l + V i3 )xH (c 3 )/i(cj) + —Q^ h ( c i) 



(21) 



(22) 



Equation (|2ip describes the one-time correlation between fluctuations in the stationary state. As it can be seen, the 
correlation function is determined by the properties of the linearized Boltzmann-Fokker-Planck operator, A, and by 
the one-particle distribution function, xh- 

For the purpose of the next section, it is also convenient to define a new function 



h H (l 12 ,Ci,c 2 ) = Xh(ci)6(1 12 )S(c 12 ) + <7if(li2,ci,c 2 ). 



(23) 



Taking into account equation (|21[) and that the term C\ 2 ■ gj^Xff(cx)<5(li 2 )#(ci 2 ) vanishes identically, the equation for 
this quantity is 



A(ci) + A(c 2 ) - C12 



/lJ?(ll2,Cl,C 2 ) 



where 



-<s(ii 2 )r( Cl ,c 2 ) + i 2 ^jf- ^ ■ ^Xh( Ci )xh(c 2 ), 



r(ci,c 2 ) =T(c 1 ,c 2 )xh(c 1 )xh(c 2 ) - [A(ci)+A(c 2 )]xh( Ci )<5(ci 2 ). 



(24) 



(25) 



C. Linearized Boltzmann-Fokker-Planck equation 



In this section we derive the evolution equation for a small perturbation around the homogeneous stationary 
distribution function (v). The objective is to evaluate some spectral properties of the linear operator that controls 
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the dynamics, which turns out to be the linearized Boltzmann-Fokker-Planck operator introduced above. These 
properties will be useful in the last section where we derive hydrodynamic equations. 
Let us introduce the small perturbation, 5f 

6f(x u t) = f 1 (x u t)-f H (v 1 ), ^-«1. (26) 

The linearized equation for 8f is obtained from the Boltzmann-Fokker-Planck equation {7J and equation (fT5|> 

Q Q 

-8f( Xu t)+v 1 ■ —Sf(x u t) = K( Xl ,t)Sf(x u t), (27) 

where K(x\,t) is defined in Q14p . Now, let us introduce a dimcnsionless perturbation and a dimensionlcss time scale 
as 



Sf{x!,t) = -^-5x(I,c,s) 

4 



(28) 



s = -t, 



which is proportional to the number of collisions per particle in the interval (0,t). The equation for 8\ reads 



^<5x(l,ci,s) = 



A ( c i) -ci • ^ 



(29) 



(30) 



where A is the linearized Boltzmann-Fokker-Planck operator defined in (|22|) . Equation (|30|) is the linearized 
Boltzmann-Fokker-Planck equation and it describes the dynamics of any small perturbation around the homoge- 
neous stationary state. 

In [32] it was shown that the linearized Boltzmann-Fokker-Planck operator has d + 1 eigenfunctions associated to 
the null eigenvalue (in principle there are not more eigenfunctions associated to this eigenvalue) . These eigenfunctions 
are 



1 d 

^( c ) = qIT ' I c Xh(c)} +xh(c) 
3 oc 



(31) 



oc 



(32) 



They fulfill A(c)£j(c) = for i = 1,2. Moreover, as the number of particles and momentum are conserved in collisions, 
we have 



J dcA(c)h(c) = J dc Cl A(c)h(c) = 0, 

where h is an arbitrary function. Equivalcntly, it can be said that the functions 

6(c) = Xff(c), £ 2 (c) = Xh{c)c, 
are left eigenfunctions of A associated to the null eigenvalue with the scalar product 

(f\g) = J dc XH \c)r(c)g(c), 

with /* the complex conjugate of /. 



(33) 



(34) 



(35) 



III. THE BOLTZMANN-LANGEVIN EQUATION 



The starting point for the study of the fluctuations in this work will be the Boltzmann-Langevin equation. In 
this section we derive the corresponding equation and we determine the properties of the noise in order to obtain 
consistency with the equations of the correlation functions presented in the previous section. 
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As defined in Eq. ([5]), the one-particle distribution function, /i is the ensemble average of the phase function Fx, 
and its dynamics is given by the Boltzmann-Fokker-Planck equation, Eq. |7|). The problem is now to find an evolution 
equation for the fluctuating quantity 

6F(l,c,s) = ^[Fx(x,t)-f H (v)}. (36) 
n H 

As for the velocity of a Brownian particle [3(| , we expect that the difference between the equation for the fluctuating 
quantity 5F, and its average, (SF) = 5x, is a random force term, R [3lj . Then, the fluctuations 5F around xh are 
described by a Boltzmann equation linearized around the solution xh with a random force, R, added 



A(0-c-| 



6F(l,c,s)+R(l,c,s). (37) 



Taking averages in equation (|37p , we obtain the linearized Boltzmann- Fokkcr-Planck equation, Eq. (|30[) , if and only 
if 

(R(l, c, s)} = 0. (38) 

Equation (|37[) is the Boltzmann-Langevin equation. As in the free cooling case [l5j . we assume that the noise term, 
R(l, c, s), is Markovian 

(R(lx,cx,sx)R{h,C2,s 2 )}x = H(h, 1 2 , ci, c 2 )S(s 1 - s 2 ), (39) 

where (. . . )h means average in the stationary homogeneous state. In order to evaluate the function H(\\, I2, Ci, c 2 ) 
cxplicitely, we will calculate (5F(lx, Ci, s)<5i ? (l 2 , c 2 , s))_sf with the Boltzmann-Langevin equation and then, we will 
impose compatibility with the equation of the correlation function of the previous section. So, let us first write this 
function as a functional of the distribution and the correlation functions. Using the definitions of the microscopic 
densities, Eqs. ([3]), Q, the distribution and correlation functions, Eqs. ([5]), (fl2)) . ([23)) . and the dimensionlcss 
distributions, Eqs. (fT6|) . (|20]l . we have 



(<5F(xi,*)<5F(x 2 ,;)) ff = -|-<[J'i(a:i,t)-/ ff (vi)][f , 1 (a; 2 ,t)-/ ff (v 2 )]) H 

„2d 



[(Fx (an, t)Fx (x a , t)) H - f H (vi)f H (v a )] 



- J t[f2,H{xi,x 2 ,t) + /ff(vi)5(a;i - x 2 ) - /h (vi)/ ff (v 2 )] 

-f-[32,i?(^i,^2) + /ff(vi)(i(a;i - x 2 )] 

[<m(li2, ci,c 2 ) + <5(li 2 )(5(ci 2 )xff(ci)] 



= ^h H (h2,cx,c 2 ). (40) 

Now, let us solve Eq. (f3"7} as a functional of the noise. In order to do that it is convenient to define the linear 
operator 

AOi.cO^AtcO-Ci-^-, (41) 
in terms of which the solution for SF(l, c, s) is 

6F(l,c,s) = e A(1 ^ c)s (5 J F 1 (l,c,0)+ [' ds'e A ^ s - s,) R(l, c, s') 

Jo 

s ^ ^ d S 'e A ^ s -^R(l,c,s'), (42) 
Jo 

where we have assumed that the term stemming from the initial conditions vanishes in the long time limit. This 
is equivalent to assuming that the spectrum of the linearized Boltzmann-Fokker-Planck operator is such that any 
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perturbation without component in the subspacc associated to the null eigenvalue decays. With equation (|42]1 . we 
can evaluate 

(5F(l 1)Cl , s)SF(h, c 2 ,s)) H 

ds' [' d s " e A(li^i)(^-^')+A(l 2 ,c 2 ) (s - s ") ( R (l uCli8 ') R (l 2i C2) S ")) H 



ds' / ds // e A(ll ' Cl)( '-*' )+A(,3 ' Ca)(a - a " ) ff(lia,ci J ca)<y(« / -s") 
Jo 

r d S ' e [ A ( ii ' ci )+ A ( i2 ' c2 )]( s - s ')i/(i 12 , Cl , c 2 ) 

Jo 



= -[A(l 1 ,c 1 )+A(l 2 ,c 2 )]- 



D [A(h,ci)+A(h.c 2 )](s-s') 



s'=0 



if(li 2 ,ci,c 2 ) 



-[A(1 X) Cl ) + A(l 2 , c 2 )]- 1 H(h 2 , ci, c 2 ), 



(43) 



where we have assumed that the term e[ A ( ll > Cl )+ A ( l2 ' C2 )] s iJ(l 12j ci, c 2 ) — > in the long time limit. After identifying the 
function, we will see that this is, in fact, the case, because iJ(li 2 , Ci, c 2 ) does not have components in the subspace 
associated to the null eigenvalue. Equivalently we have 



[A(li,ci) + A(l 2 ,c 2 )]^(li 2 ,ci,c 2 ) = -n^ d ir(li 2) ci,c 2 ). 
Finally, comparing equations (|4"4"|) and (|24[) we conclude that 

1 | 2 d d 

i?(li 2 ,ci,c 2 ) = — ^5(li 2 )r(ci,c 2 ) - j^q^ ■ q^-Xh(ci)xh{c2), 



(44) 



(45) 



with r(ci,c 2 ) given in Eq. (|25|) . With this expression of H we can see that it does not have components in the 
subspace associated to the null eigenvalue. Taking into account that £i and £ 2 are left eigenfunctions of A associated 
to the null eigenvalue and that [32| 



dci J dc 2 T(ci,c 2 )xij(ci)xff(c 2 ) = 0, 
dci I dc 2 ciiC 2i T(ci,c 2 )xif(ci)xH(c 2 ) = | 2 , 



(46) 



we can see that 



J di i2 J dc x J dc 2 -ff(li 2 ,ci,c 2 ) = 0, 
J d\\2 J dci J dc 2 ciiC 2i iJ(li 2 ,Ci,c 2 ) = 



(47) 



In the remaining of this section, we shall write the Boltzmann-Langevin equation together with noise properties 
in Fourier space. This will prove useful for the subsequent analysis. Let us introduce the Fourier component of a 
function of the position variable as 



/(k) = J die-* 1 / (1), /(I) = i e<k l /»' 



where V = jj is the volume in units of the mean free path. The equation for 5F(k, c, s) is then 



d_ 

ds 



A(k,c) 



<5F(k,c, s) = ^(k,c, s), 



where we have introduced the operator 

A(k, c) = A(c) - ik ■ c. 
The Fourier transform of the noise, i?(k, c, s), obeys 

(R(k,c,s)) H = 0, 



(48) 

(49) 

(50) 
(51) 
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and 



where 



(R(ki,a, s 1 )R(k.2,c 2 , s 2 )}h = H(k 1 ,k 2 ,c 1 ,c 2 )S(s 1 - s 2 ), 



i7(ki,k 2 ,ci,c 2 ) 



N 



r(ci,c 2 )J(ki +k 2 ) -C 2 ^ • ^Xff(ci)xff(c 2 )<5(ki)5(k 2 ; 



(52) 



(53) 



which completes the calculation of the noise variance within the Langevin description, that will be used to quantify 
the transverse velocity fluctuations. 



IV. THE FLUCTUATING TRANSVERSE VELOCITY 



The objective in this section is to derive a fluctuating equation for the transverse velocity field, Wj_(k, s). The 
reason to consider this field is that its equation is decoupled from the equations for the other hydrodynamic fields 
[2(i| , and it can be derived exactly in the hydrodynamic limit with the knowledge we have about the spectrum of the 
linearized Boltzmann-Fokker-Planck operator. 

Mathematically, the transverse velocity is defined in the following way: Let us consider the d-dimensional vector k 

d-l 



which belongs to the space 5R d . This space can be expanded in the base |k| (J |k^| where k is a unitary vector 
parallel to k and |kj_ j are (d — 1) unitary vectors orthogonal to k. The transverse velocity is defined as 



w±i(k, s) — / dc(c ■ k l ± )SF(k, c, s), 



,,d- 1. 



For the subsequent analysis, is is convenient to introduce the following projectors 

P»/(c) = <6±i(c)|/(c))^(c), 



)<<> 



/(c) = (1-PW)/(c) 



Then, if we apply PW to SF we obtain 



P^6F(k,c,s)=w ±i (k,s^ 2±i (c), 



(54) 

(55) 
(56) 

(57) 



and the transverse velocity is the component of <5F(k, c, s) in the subspace generated by £ 2 ±;(c). 

In order to obtain an equation for wi, we apply the projectors P and Q (for simplicity in the notation we skip the 
super-index) to the Langevin equation P^)) 



|"-PA(k,c) 

OS 



PSF(k, c, s) = PR(k, c, s) - Pi(k ■ c)Q6F(k, c, s), 



(58) 



— -QA(k,c) 

OS 



Q6F(k, c, s) = QR(k, c, s) - Qi(k ■ c)P5F(k, c, s) 



(59) 



Now, let us solve equation (|59")) formally 

Q5F(k,c,s) = e QA( - k ^ s Q5F(k,c,0) 

+ f ds'e QHKc)(s ~ s,) [QR(k,c,s')-Qi{k-c)PSF(k,c,s')}. 
Jq 

In the time regime in which the system has forgotten the initial condition, i.e. when we can consider that 
e QA( - k -'°^ s QSF(k, c, 0) -> 0, by substituting Eq. (60]) in Eq. J5HJ), we obtain a closed equation for PSF(k, c, s) 



(60) 



|--PA(k,c) 



P5F(k,c,s) + P(k-c) f ds'e QA( - k - c ^ s - s ">Q(k-c)P5F{k,c,s') 
Jo 

= PR(k,c,s) 

- Pi(k-c) I ds'e QA{k ' c){s - s ' } QR(k,c,s'). 
Jo 



(61) 
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In Appendix A it is shown that, in the hydrodynamic limit, i.e. to second order in k and in the long time limit, 
equation (|61[) reduces to the following equation for the transverse velocity 



os 



w±(k, s) = R w (k, s) 



The coefficient fj is the shear viscosity given by 



n 



J ds J dcc x c y e A( - c)s c x £ 2 , y (c) = - J dcc x c y A(c) l c x £ 2 ,y(c), 



(62) 



(63) 



which agrees with the expression obtained in [2J| by the Chapman-Enskog method, and R w (k, s) is a noise term which 
can be decomposed as 



i4,(k,s) = 5(k,s) + n(k, s). 



(64) 



The term SQs.,s) comes from the thermostat (which does not conserve momentum locally), and the second term, 
II(k, s), is the fluctuating part of the pressure tensor. Their microscopic expressions in terms of the noise of the 
Boltzmann-Langevin equation are 



and 



5(k,s) = / dc(k ± ■ c)R(k,c,s), 



QA(k,c)(s-a') 



n(k, s) = -ik / ds' / dc(k ■ c)(k ± ■ c)e 



QR(k,c,s'). 

Of course, due to the fact that (R(k, c, s))h = 0, the mean value of the noise vanishes 

(R w (Ks))h = 0. 



(65) 



(66) 



(67) 



The autocorrelation function of the noise is evaluated in Appendix B. Due to symmetry considerations, there are only 
correlations between the k and — k components, which yields 



(Rw(k, si)R w {-k, s 2 ))h 



N 



2 . 



i 2 6(si - s 2 ) + k 2 C xy (s 2 - si) 



(5(k, Sl )n(-k, s 2 )) H , 



Sl < s 2 . 



(68) 



The first term is the expected zeroth order term which comes from the heating. The Dirac delta function is an exact 
consequence of the fact that the external noise (the heating) is white. The function C xy (s 2 — s±) reads 



C xy (s 2 — si) 



dc! I dc 2 c lx c ly c 2x c 2 ye h{c2){s2 Sl) cj) H (c u c 2 ). 



(69) 



Here, we have introduced the function 4>h(ci,c 2 ) as the space integral of the correlation function hn (1, ci, C2) 

0ff(ci,c 2 ) = J d\h H {\, ci,c 2 ) = xh(ci)5(ci 2 ) +X2(ci,c 2 ), (70) 

where X2( c ij c 2 ) = / ^15^(1, Ci, c 2 ) and gn is the dimensionless correlation function defined in (|20p . The equation for 
can easily be obtained by integration over space variable in Eq. (|24[) . which gives 

(71) 



[A(ci)+A(c 2 )]^(ci,c 2 ) = -r( Cl ,c 2 ) + e 2 ^ • q^xh(ci)xh(c 2 ). 



As discussed in Appendix B, C xy can be physically interpreted as the autocorrelation function of the global quantity 
Yli—i Vx{t)Vy{t)- in the elastic case, this correlation function is related to the shear viscosity but it is not the case for 
granular systems [Hj]. The formula for the correlation (5(k, si)IT(— k, s 2 ))h is given in Appendix B, and it does not 
seem to admit a direct physical interpretation. As the two correlation functions, C xy {s) and (S(k, Si)II(— k, s 2 ))fl-, 
decay with the kinetic modes and, in the k — > limit, the fluctuating velocity is expected to be frozen (its time 
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evolution is given by the null eigenvalue, Eq. (J62J) , we can consider that they are proportional to a Dirac delta 
function in time, and we have 



^-k 2 C xy (s 2 - si) + (5(k, fli)n(-k, s 2 )) H 
~V 2 

—k 2 

N 



dsC xy {s)+ / <is(S(k,0)II(-k,s))ij 
Jo 



S{si - s 2 ). 



(72) 



Note that this is in contrast with the free cooling case where the noise can not be considered to be white [l5j. In 
this case, the equation for the transverse velocity (rescaled with the thermal velocity) contains a term of order zero 
in k which is proportional to the cooling rate. Then, even in the k — > limit, the velocity is not frozen on the kinetic 
scale. In Appendix C the second integral of equation (|72|) is evaluated obtaining 



(73) 



J ds(S(k, 0)II(-k, s)) H -> ~Jj~ k2 j dCl j dc 2Cl x C 2x C 2y A(c 2 ) 1 c 2y (j)H (ci,c 2 ), 



which is valid in the hydrodynamic limit. If we substitute 4>h{ci, c 2) by its expression in terms of the one and two- 
particle distribution function, Eq. (|70[) . wc obtain (see Appendix D) that the one-particle contribution vanishes and 
the correlation function can be written in terms of the two-particle velocity correlation function, X2(ci, C2) 

(R w (k, si)Ru,(-k, s 2 ))h 



^S( Sl -s 2 )< 



- J dci J dc 2 c lx c ly c 2x c 2y A(c 2 ) ^(ci, c 2 ) 
j dct J dc 2 c lx c 2x c 2y A(c 2 y 1 c 2y x 2 {c 1 ,c 2 ) j . 



(74) 



As i? u ,(k,Sx) is Gaussian, the noise is completely characterized by Equations (|6"T|) and (jT4"|) . As it can be seen, the 
k 2 term has no relation, a priori, with the shear viscosity, Eq. (|63p . i.e. there is a priori no fluctuation-dissipation 
relation as that assumed in (20| . However, we now show that, under additional hypothesis (that in principle are not 
restricted to the elastic limit), the aforementioned term reduces to the shear viscosity. Let us assume that the most 
important contribution of the two particle velocity correlation function, x 2 (ci,c 2 ), is the hydrodynamic part, i.e. we 



d+2 d+2 

X 2 (C1,C 2 ) ~ ^ «/3:/3'£/3( C l)£/3'( C 2)- 
/3=1/3'=1 



(75) 



This assumption was already made in [32| where the coefficients ap t p> were evaluated to calculate the total energy 
fluctuations. We emphasize that it led to an excellent agreement between analytical predictions and numerical data 
(Monte Carlo) for energy fluctuations, for all the values of the inelasticity [32 |. In this approximation, the first integral 
in (|74p vanishes because c x c v xh (c) is orthogonal to the hydrodynamic modes. The second term is 



/ dci / dc 2 ci x c 2x c 2y A(c 2 ) 1 c 2;/ y^ ^ ap,p'£p(ci)£/3'(c 2 ) 
J J p=i/3'=i 

= a 2xax J dc lC ^ 2x ( Cl )J dc 2 c 2x c 2 yA(c 2 )-^ 2x (c 2 ), 



(76) 



where, for symmetry considerations, the only term that remains is the one associated to (3 = (3' = 2. If we use now 
that a 2 i 2 i = — 1/2 (see the reference [HJ), and the formula for the shear viscosity, Eq. (|63|) . we finally have 



(R w (k, si)i? lu (-k, s 2 ) 



N 



5{si-s 2 )(¥+r}k') 



(77) 



Then, in the hydrodynamic limit and assuming that the two-particle velocity correlation function, X2(ci,c 2 ), has 
only components in the hydrodynamic subspace, the correlation function of the noise reduces to the one introduced 
phcnomcnologically in 20j, not only in the clastic limit, but for arbitrary inelasticity Although we have no direct 
proof of the accuracy of approximation (|75p . wc note that it is backed up by numerical data, see e.g. [32| . 
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V. CONCLUSIONS AND OUTLOOK 



The primary objective of this work was to derive a Boltzmann-Langevin description for a heated granular gas, in 
the spirit of the work of Bixon and Zwanzig [3l| . who considered conservative systems. The system considered here 
is on the other hand dissipative, and is heated by a stochastic force changing the velocity of the particles between 
collisions. The loss of energy in collisions is compensated by the energy given to the particles by the thermostat and 
a stationary state is thereby reached. Our study is restricted to this stationary state. The Boltzmann- Langevin 
equation has been derived and the properties of the noise appearing in this equation were identified by imposing 
consistency with the equation for the correlation function derived in [32j | . 

The Boltzmann-Langevin equation is the starting point for the derivation of fluctuating hydrodynamic equations. 
This can be done by projecting the Boltzmann-Langevin equation into the hydrodynamic subspace. As our knowledge 
of the spectrum of the linearized Boltzmann- Fokker-Planck operator is quite limited, we have focused on the equation 
for the transverse velocity field, that is decoupled from the rest of the fluctuating hydrodynamic equations. This 
specific case was studied in [Til [l5j for the free cooling state, where it was shown that the relevant Langevin noise 
is not white and that there is no a fluctuation-dissipation relation. In other words, the amplitude of the noise is not 
related to the shear viscosity. On the other hand, in the stochastically heated system, the behavior that we have 
reported is different. First, the noise of the transverse velocity contains two parts: one coming from the thermostat 
(which does not conserve momentum locally) together with the more standard fluctuating pressure tensor. The 
correlation function of the noise can then be written as a sum of direct and cross terms made up from the previous 
two contributions. Moreover, in contrast to the free cooling scenario, the noise can be considered as white J3f|, as 
the dynamics of the velocity is as slow as desired in the hydrodynamic (low k) limit. As in the free cooling case, in 
the hydrodynamic limit, the amplitude of the noise is priori not related to the shear viscosity (such a relation, of 
fluctuation-dissipation type, had been assumed in the approach of Ref [I(|). However, considering that the two-particle 
velocity correlation function has only hydrodynamic modes -which seems a reasonable assumption- somehow restores 
fluctuation-dissipation and we obtain the expression assumed in [2fj| . with the actual inelastic shear viscosity. In 
principle, this rather surprising result - reminiscent of those reported in Refs (37j- is not limited to small inelasticity. 

For future work, remains the extension of the theory to the other hydrodynamic equations (beyond the transverse 
velocity), together with a generalization of the scheme presented here to a more general class of thermostated systems, 
such as the ones with multiplicative noise [35| . 
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Appendix A: DERIVATION OF THE TRANSVERSE VELOCITY FIELD EQUATION 



In this Appendix we derive the equation for the transverse velocity field, u>j_(k, s), in the hydrodynamic limit. The 
starting point is equation (|6Tj) 



£-PA(k,c) 



PSF(k,c,s) + P(k-c) f ds'e QA ^ c K s - s ^Q(k-c)P5F(k,c,s') 
Jo 

= PR(k,c,s) 

- Pi(k-c) [ S ds'e QA( - k ' c ^ s - s,) QR(k,c,s'). 
Jo 



(Al) 



Let us first consider the term PA(k, c)P8F(k, c, s). As A(c)£2_i_(c) = and J dc(k± ■ c) (k • c)£2_L (c) = 0, we easily 
have 



PA(k, c)PSF(k, c, s) = 0. 



(A2) 
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Let us evaluate the last term of the left hand side of Eq. (|A1|) . To second order in k, we have 

P(k-c) f ds'e QA( - k ' c ^ s - s '^Q(k-c)P6F(\s.,c,s') 
Jo 

^fc 2 6±(c)^ ds'wx(k,s') y*dc(k-c)(k ± •c)e A(c)( ^ s ' ) k-c6±(c) 

= fc 2 6±(c) I" ds'w x (Ks')G xy (s-s') (A3) 
Jo 



where we have introduced 



G xy ( S )= /" dc(k ± . C )(k-c)e A ( c ) s (k-c)6 ± (c) 

dc C;cSe A ( c ) s c ^2,(c), (A4) 



and use has been made of the fact that the operator A(c) is isotropic. In the hydrodynamic limit, the velocity evolves 
in a scale much slower that the scale in which the function G xy (s) decays. We then have 



ds'w±(k, s')G xy (s - s') -* fjw±(k, s), (A5) 
where fj is the dimensionless shear viscosity 

fj= / dsG xy {s). (A6) 



o 



The noise terms are the last two terms of Eq. (|A1|) 

PR(k,c,s) = 6±(c) J dc(k x -c)R(k,c,s) 



= 6±(c)5(k, s ), (A7) 

and 

P(ik ■ c) f ds'e QA{Kc)is - s,) QR(k, c, s') 
Jo 

= &x(c)ik J dc(k_L-c)(k-c) / ds'e QA( - k ' c ^ s - s '^QR(k,c,s') 

= -£ 2 ±(c)II(k, s ), (A8) 

where we have used the definitions of S(k, s) and the fluctuating pressure tensor, II(k, s), Eqs. ([65]) and (|66|) . Finally, 
by multiplying Eq. (|A1[) by k± ■ c and further integrating over velocities, we obtain the equation of the transverse 
velocity of the main text. 

Appendix B: AUTOCORRELATION FUNCTION OF R w (k,s) 

In this Appendix we evaluate the correlation function of the noise of the transverse velocity field, R w (k,s). We 
consider k ^ 0, si < S2 with si large. It is convenient to introduce the following notation for the transverse and 
parallel components of the vector c 

k± ■ c = c±, k-c = C||. (Bl) 

The autocorrelation function of R w (k, s) reads, in terms of S(k, s) and II(k, s), 

(R w (k, Sl )R w (-K s 2 ))h = {(S(k, Sl ) + n(k, «i))(5(-k, s 2 ) + n(-k, s 2 )))h 

= (S(k, si)5(-k, 8 2 ))h + (S(k, Sl )n(-k, s 2 )) H 

+ (n(k, Sl )S(-k, aa)) h + <n(k, si)n(-k, S2 )) H . 

(B2) 
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N 

-2 



We now calculate each correlation function taking into account the microscopic expressions of S(k, s) and II(k, s), 
Eqs. (frJ5)) and ([fjFJ|) , and the correlation function of the noise of the Boltzmann-Langevin equation, Eqs. and (|53p . 
The first term is 

(S(k,s 1 )S(-k,s 2 ))H = J dc! J dc 2 c 1 ±c 2± (R(k,c 1 ,s 1 )R(-k,c 2 ,s 2 )) H 

S(si - s 2 ) J da J dc 2 ci_LC2_i_r(ci, c 2 ) 
— 5(si-s 2 ) J dci /'dc 2 ci_ L c 2 _LT(ci,c 2 )x J j(ci)xjj(c 2 ) 

= f ^( fll - aa ). (B3) 

where we have used the relation J dci J dc 2 ci±c 2 ±T(ci,c 2 )xH{ci)xH{c 2 ) = £ 2 , that is proved in [32| . 
The second correlation function is 

(S(k, Sl )U(-k,s 2 )} H 

= J dc lCl± ik J ' ds J dc 2 c 2ll C2±e Q ^- k ' c ^-^(R(k,c 1 ,s 1 )Q 2 R(-k,c 2 ,s)}H 

= lk Yl J dci y"dc 2 c 1± c 2 || C2± e Q2A (- k ^» S2 - Sl )Q 2 r(c 1 ,c 2 ). (B4) 

Here, we have changed the sign of II(— k, s) because we are dealing with the — k component and then k — > — k (we do 
not change kx — > — k±, because this vector comes from the projector P and it is fixed). 
The third term vanishes 

(II(k, si)S(— k, s 2 ))h = —ik J ds' J dc 1 c 1 \\a± J dc 2 c 2 ±_ 

x e QlA( - k ' c ^ s -^Q 1 (R(k,c 1 ,s')R(-k,c 2 ,s 2 )} H = 0, 



(B5) 



because (R(s')R(s 2 ))jj = for s' € (0,si) with Si < s 2 . 
Finally, we evaluate the last term to second order in k 



(n(k, Sl )n(-k,s 2 )> H 

~fc 2 jT ds[ j S \s' 2 J d Cl J dc 2Cl || Cl±C2 ||C 2±e A ( Cl » Sl - s ^ +A ( C2 )^- s ^ 



x(i?(k, Cl ,si)i?(-k, c 2 ,s' 2 )) H 
>2 r s i 



= Yfk 2 J* ds[ jj ds' 2 J dci y ^Cijic^c.jic^e^ 1 )^-^)^)^-^^^,^)^ -4) 
= ^-k 2 J d Cl y dca^iicixcaiicaxy^ds^J^-'H^J^-'JrCci.ca) 
= -^fc 2 | dd y rfc 2Cl || Cl± c 2 || C2i e A ( Cl ^ +A ( C2 ) S2 £ d S e- s ^+ A ^r( Cl ,c 2 ) 



—k 2 I d Cl / dc 2Cl || Cl±C2 ,|C 2±e A ( C5 



TV 

x[A( Cl )+A(c 2 )]- : 



- d d 
-r(c l5 c 2 ) + i 2 — ■ ^—Xh(ci)xh(c2) 
oci ac 2 



(B6) 



where, in the last step, we have taken into account that si is large and we have introduced the term ■ 
-^Xh(ci)xh (c 2 ). This term does not contribute to the integral, but is written for convenience, to make a con- 
nection with the global correlation function (j>H (ci, c 2 ) = J dl/i#(l, ci, c 2 ) which fulfills Eq. (fTTjl . In doing so, we find 
that the autocorrelation function of II(k, s) reads 

(II(k, Sl )n(-k, s 2 )) H ~ ^r^C^Csa - si), (B7) 
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where 



C xy (s 2 - si) = J dci J dc 2 c lx c ly C2xC2ye A2( - S2 Sl) (t)H(ci,c 2 ). (B8) 
As 0ij(ci, C2) is the integral of the correlation function hn (1, ci, C2), can be identified as the correlation function 



of the global quantity J2i=i ^(t)V y (t) 



Appendix C 



In this Appendix we evaluate the time integral of the correlation function (5(k, 0)II(— k, s))h i n the hydrodynamic 
limit. Using the notation of the previous Appendix, we have 



ds(5(k,0)n(-k,s)) Jf 







V 



ik— I dci I dc 2 c 1 ±c 2 \\C2 



= Q 2 A(-k,c 2 )s 



Q 2 r(ci,c 2 ) 



> 2 A(-k,c 2 
V 2 f f 1 

ik ^7 I rfc l / dc 2 C 1± C 2 uC 2± — rQ2r(ci,C 2 ) 

N J J Q 2 A(-k,c 2 ) 





2 



Q 2 A(-k,c 2 ) 

.,V 2 f f Q 2 [A(k, Cl ) + A(-k,c 2 )] ^ 

ik— / dci / dc 2 c 1± c n c 2± r </> ff (k,Ci,c 2 ), 

A 7 7 <y 2 A(-k,c 2 ) 



(CI) 



where we have introduced the function 



0H(k,ci jCa ) = J d\e~^ l h H (\,c l ,c 2 ), (C2) 
that fulfills the Fourier transform of Eq. ([24]) 



:2 a 9 



[A(k,ci)+A(-k I c a )]0 jr (k > ci,ca) = -r( Cl ,c 2 )+^— ■ t— Xh (ci) X h (<3>)*(k). (C3) 

aci ac 2 

Note that the last term in the previous equation only appears for k = 0. Taking into account that c±\h{c) is left 
cigenfunction associated to the null eigenvalue and after some algebra we obtain 

ds(S(k,0)U(-k,s)) H 

-k 2 [ dc! I dc 2 c 1 nc 1± c 2 nc 2± 1 r0 ff (k,Ci,C 2 ) 

T J J Q 2 A(-k,c 2 ) 

ik J dci J rfc 2 c 1 _ L c 2 ||C2_L0fl-(k, ci,c 2 ). (C4) 



N 

Now, let us consider the hydrodynamic limit of (|C4[) . The first term gives 



V 2 f f 1 V" 2 f°° 

— k 2 J dci y dc2Ci||Ci_ L c 2 ||C2_L^-^-j^-y</)ff(k,c 1 ,C2) -> - — fc 2 y dsC xy (s). (C5) 



The second term can be evaluated by using the following expansion in powers of k 

[A(c) - ik ■ of 1 A(c)" 1 + A^)" 1 ^ • c)A(c)" 1 , (C6) 

which yields 

U 2 9 r 00 , » v 2 



N 



ik J dci y dc 2 ci ± c 2 ||C2±0_ff(k, ci,c 2 ) 

f°° V 2 f f f°° 

k 2 J dsC xy (s) ~ —k 2 J dcx J dc 2 c lx c 2x C2y J dse Mc2)s c 2y (j) H (ci, c 2 ). 



(C7) 
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Taking into account (|C5|) and (|C7[) . wc obtain 

>2 



J ds(S(k,0)n(-k,s)) -> -— fc 2 y dcx y dc 2 c lx c 2x c 2y J dse K{ - C2)s c 2y (j) H {c u c 2 ) 

= -^-fc 2 y rfci J dc 2 c lx c 2x c 2y A(c 2 y 1 c 2y (f>H{ci 7 c 2 ). (C8) 

Appendix D 

In this Appendix we evaluate the fc 2 component of the correlation function of -Ru, in terms of the one and two-particle 
distribution functions, xh and X2 

<Aff(ci,c 2 ) = Xh(ci)5(ci 2 ) +X2(ci,c 2 ). (Dl) 

The first term is 

y dsC xy (s) = J ds J dci J dc 2 ci x ciyC 2x c 2 ye s ^ C2) <pH{ci,c 2 ) 
= - j dci J dc 2 ci x ci y c 2x c 2y A(c 2 )~ 1 (t>H(ci,c 2 ) 
= - J dccxCy^cj^CxCyXHic) 

- j ' dci J dc 2 ci x ci y c 2x c 2y A(c 2 )~ 1 X2(ci,c 2 ). (D2) 
If we do the same in the second term, we have 

J dci J dc 2 c lx c 2x c 2y A(c 2 y 1 c 2y (l) H (c 1 ,c 2 ) = J dcc x c y A(c)~ 1 c x c y XH (c) 

+ J dci J dc 2 c lx c 2x c 2y A(c 2 y 1 c 2y X2{ci,c 2 ). 

(D3) 

It can be seen that the sum of the two terms only depends on the two-particle correlation function and wc obtain the 
k 2 part of Eq. (J71J). 
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